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We consider the emission spectrum of high-energy electrons in an intense laser field. At high 
intensities (ao ~ 200) we find that the QED theory predicts a narrower angular spread of emissions 
than the classical theory. This is due to the classical theory overestimating the energy loss of the 
particles, resulting in them becoming more susceptible to reflection in the laser pulse. 


I. INTRODUCTION 

Since the discovery of chirped pulse amplification [T] 
the powers and intensities of state of the art laser facil¬ 
ities have been exponentially increasing [2], the current 
record of 2 x 10 22 W cm -2 having been set in 2008 [3]. 
With the advent of various new facilities over the next 
few years, such as the Vulcan 20 PW upgrade [4j, the 
Extreme Light Infrastructure (ELI) Facility [5] and the 
XCELS project |6], this trend is expected to continue 
for the foreseeable future. The availability of such tech¬ 
nology has driven a large field of research in the topic 
of nonlinear Thomson and Compton scattering, the un¬ 
derstanding of which is important from a fundamental 
physics perspective 130 . Additionally, and at least as 
importantly, the process produces high-energy, tuneable 
y-ray beams, which are important for fundamental re¬ 
search [9], as well as more practical applications such 
as cancer radiotherapy m and the radiography of dense 
objects El- Recent experiments P2UE9 have been push¬ 
ing the limits of peak energies and brilliances, taking us 
towards the regime where radiation reaction and QED 
effects will start to come into play [141116] . (For related 
studies of high energy electrons in orientated crystals see, 
e.g., [TTHU] ) 

In this article we study nonlinear Thomson and Comp¬ 
ton scattering at ultra high intensities, assessing the im¬ 
pact of classical radiation reaction and QED effects on 
the properties of the emitted photon spectra. We find 
that, because the classical theory overestimates the ra¬ 
diative energy loss, it predicts a broader angular spread 
of emissions than the QED theory. 


II. THEORY 

We consider the case of an electron in a head on col¬ 
lision with a laser pulse. We adopt natural units where 
h = c = 1. To begin with we will take our laser to 
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be a plane wave field propagating in the ^-direction de¬ 
scribed by the null wave vector k M = u;o(l,0, 0,1), with 
central frequency u)q. We assume the field to be po¬ 
larised in the perpendicular (x) direction and therefore 
introduce the polarisation vector e = (0,1,0,0). We 
define the dimensionless intensity in the usual manner, 
ao = eE/uj^m, where E is the peak magnitude of the 
electrical field strength. The electromagnetic field ten¬ 
sor of the wave is taken to depend arbitrarily on the 
phase (f> = k • x = uo{t — z); F^ u ((j)) = ao/(0)/ Mzy , where 
f^ v — (k^e 1 ' — k v e^)/uj and /(</>) is a function describing 
the pulse. Additionally, to aid future discussion, we also 
define the time-dependent intensity a(</>) = ef((j))E/ujm , 
which is equal to ao at the pulse peak. 

The emission spectrum from a particle in the field can 
be decomposed into a sum of harmonics, corresponding 
to multiples of the laser frequency. In the quantum de¬ 
scription these correspond to the number of laser photons 
involved in the scattering process m- During each scat¬ 
tering process, the electron will absorb an integer number 
n of laser photons (each of momentum k) before emitting 
a photon of momentum k' . From conservation of momen¬ 
tum arguments it can be shown that the frequency of the 
scattered photon is given by [21 J 


where 


nujQ 

1 + jn(l - COS#) ’ 


nujo/m - 7/3 + a§7(l - /?)/2 

7(1 + £) 


(1) 

( 2 ) 


It can be seen that when j n < 0 the maximum emis¬ 
sion frequency occurs when the photons are backscat- 
tered (0 = 180°). Conversely, when j n > 0 the max¬ 
imum frequency occurs for forward scattering (0 = 0). 
Note that the support for a given harmonic depends on 
the harmonic number n (and that the amplitude of a 
given harmonic will be determined by the cross sections 
provided in the references). For high intensities the spec¬ 
trum will be composed of a very large number of harmon¬ 
ics, with the spectrum decaying after the harmonic with 
number n ~ 3/2 [22J. 

Now let us consider the angular directions of the emis¬ 
sions. From expressions 0 and © it can be seen that 
the angular range of each harmonic will depend on both 
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FIG. 1 . Diagrams showing the relationship between the laser 
momentum and the peak scattering angle of the emitted pho¬ 
tons. (a) ao <C 7, the peak emissions will be forward scattered 
relative to the laser axis, (b) ao ~ 27 defines the ‘center-of- 
mass’ frame for the collision where the peak emissions will be 
at 90 °. (c) ao 7 the peak emissions will be back scattered 
relative to the laser axis. See Ref. m for further details. 


ao and the 7 -factor of the electron. It can be shown that 
the angle at which the emissions peak depends on the 
ratio of the sum of the laser photon momenta to the elec¬ 
tron momenta (see m for further details). In the case 
where the electron momentum is greater than the sum of 
the laser photon momenta for all the photons involved in 
even the highest order scattering processes, i.e. ao 7 , 
the peak emissions will be forward scattered relative to 
the laser axis. In the opposite case, ao ^ 7 they will be 
backscattered, and the case ao ^ 27 defines the ‘center- 
of-mass’ frame for the collision (see EUI23]) where the 
peak emissions will be at 90°. This is demonstrated in 
Fig. [TJ 

When the background field is of high intensity the ra¬ 
diation emissions from the particle will be so strong that 
the resulting energy loss will begin to effect its motion 
[22]. This radiation reaction (RR) effect will cause the 
7 -factor of the particle to decrease during the interaction 
with the laser field. This will mean that in order to esti¬ 
mate the peak emission angle we must consider the value 
of the 7 -factor when the particle is in the peak of the 
laser field (since the highest intensity region will produce 
the strongest emission signal strength and therefore dom¬ 
inate the total spectrum), rather than the initial value. 
This is illustrated in Fig. [2] where we show how the ratio 
of laser intensity to particle 7 -factor changes with time. 

If we wish to consider quantum effects then it is in¬ 
structive to introduce the dimensionles s and in variant 
‘quantum efficiency’ parameter y e = yj'p^F^ u p u / m 2 ~ 

^E/E cr , where E cv = 1.3 x 10 16 Vcm -1 is the QED ‘crit¬ 
ical’ field (‘Sauter-Schwinger’ field) [25]. This can be in¬ 
terpreted as the work done on the electron by the laser 
field over the distance of a Compton wavelength. In the 
regime y e > 1 quantum effects, including pair produc¬ 
tion from the emitted photons, will dominate. For the 
purposes of this study we restrict ourselves to the regime 
where ao, 7 1 , such that quantum effects play a role 

in the Compton spectra, but have y e < 1, so that pair 
production can be neglected. 



FIG. 2. Diagram showing how the 7-factor of the particle 
decreases due to RR loses throughout its interaction with the 
laser pulse. The dynamics can be roughly divided into three 
regimes: (i) Initially 7 a(<f>), (ii) as the particle approaches 

the peak field it loses energy due to RR, resulting in 7 ~ a( 0 ), 
(iii) after the peak field there will be a period where 7 < a(<f>). 


III. NUMERICAL TECHNIQUES 

To generate the emission spectra we use the single par¬ 
ticle code SIMLA [26], which runs in both classical and 
QED modes. In the classical case we propagate the par¬ 
ticle through the fields using the Landau Lifshitz (LL) 
equation, which takes into accounts RR effects via the 
inclusion of some correctional terms to the Lorentz force 
equation m 

=F^u„ + r 0 (F fJ - , 'u v + 

F^F a v u v - u a F av F v % , { °^ 

where ro = ( 2 / 3 )e 2 / 47 rra is the classical electron radius, 
and the overdot denotes differentiation with respect to 
proper time. We note that there are numerous alternative 
equations in the literature (for an overview see Refs. [28, 
129] and for further discussions of the LL equation see 
Refs. [30/31]). However, the LL equation has, along with 
some others, recently been shown to be consistent with 
QED to the order of the fine-structure constant a m- 

The resulting classical emission spectra are calculated 
using a novel Monte Carlo method introduced in [33]. 
The method is simple and computationally efficient. As 
the particles in our simulation are ultra-relativistic the 
radiation due to transverse acceleration is dominant, 
since this is a factor y 2 larger than that due to longitudi¬ 
nal acceleration [34] . Since the acceleration and velocity 
of the particle are perpendicular, we can approximate 
the radiation as synchrotron radiation. In our method 
we calculate the effective magnetic field , acting on 

the particle over each timestep. This is the magnetic field 
which would cause the same acceleration as the electric 
and magnetic fields together. The typical frequency of 
synchrotron emission is then given by uj c = 3eH e ^ 2 /2m. 
For an ultra-relativistic particle in an external, homoge¬ 
nous magnetic field, the classical radiation cross section 














3 


can be expressed in terms of the intensity given by [3?| 


did _ 1 31 

duo' l o' du r 


e 3 i7 eff 
27r uj'm 




(4) 


where F\(£) = £ /^°° i^ 5 / 3 (£ / )d£ / is the first synchrotron 
function. (We note that 0 is integrable in the limit u/ —>• 
0 and therefore the expression is well-defined. For further 
details see Ref. m-) At each timestep the quantity uj c is 
calculated and a Monte-Carlo method is used to sample 
from the spectra. The direction of the emission is taken 
to be that of the particle velocity, a good approximation 
for the ultra-relativistic case [34]. 

In the QED case our code works as follows. We begin 
by introducing the (quantum) synchrotron parameteiQin 
terms of the frequency of the emitted radiation a/ 


S — o / * 

6Xe — W 

In the case of high energy emissions, £ >> 1, it can be 
shown that the photon coherence length goes like [36] 

(6) 

where A = 2tt/ujq is the laser wavelength. Thus in the 
high-intensity limit ao 1 the size of the radiation 
formation region is much smaller than the laser wave¬ 
length cm and so the laser background field can be ap¬ 
proximated as locally constant and crossed [37|. We are 
then able to determine the probability of photon emis¬ 
sion using the expression for the constant crossed field 
rate T q per unit time, 


dTq 

dx 7 


V^TT-fXe 


2 + 


1 + X 


POO 

^ 2 / 3 (x)-/ d yK 1/3 (y) 

J Y 


(7) 


where K v is the modified Bessel function, x 7 = 
yJk'^F^ u k fu /m 1 2 for the emitted photon with momentum 

note that x = x 7 /(Xe ~ X 7 )> and X = 2x/( 3 Xe)- 
Although dT q /dx 7 diverges at small y 7 , the total rate 
of photon emission T q , given by integrating ^ over all 
X 7 £ [0, Xe] ^ is finite. (This apparent softening of the 
usual infra-red divergence in QED is explained in [38].) 
A discussion of the relationship between quantum expres¬ 
sion for the emission rate ^ and the classical expression 
0 is given in the appendix. 

In the QED simulations the electron is propagated 
along a classical trajectory divided into discrete time 


1 Note that in the classical regime the energy of the emitted radi¬ 
ation is typically much lower than the electron energy my (since 

the radiation is emitted continuously rather than as discrete 
emissions) and so we can approximate £ « (2/3xe)(k//7)- Fur¬ 
ther approximating Xe ~ 7 H e ^e/m 2 , one can see that £ a//a; c , 

the argument of our classical expression 



FIG. 3. Integrated frequency spectra (left panel) and angular 
spectra (right panel) for the case of an electron with initial 
70 = 800 colliding with a laser of intensity ao = 200, A = 
0.8/im and duration 30fs. Blue line: classical. Red line: QED. 


steps. After each step At the code calls a statistical sub¬ 
routine to calculate the probability of photon emission 
and to correct the electron’s momentum. The routine 
generates a uniform random number r G [0,1] and, if the 
condition r < T q At is satisfied, an emission is deemed 
to occur (under the requirement r q At < 1 ). Note that 
dr q /dy 7 (and T q ) are time-dependent quantities, due to 
the temporal variation of both the laser pulse and the 
electron motion. If an emission event occurs, a second 
uniform random number £ G [0,1] is generated and the 
photon’s x 7 (and hence its frequency) is calculated as the 
root of the sampling equation 



The photon momentum is then determined by x 7 to¬ 
gether with the assumption that at high 7 the photon 
is emitted in the direction of the electron’s motion, just 
as in the classical method described above. Finally, the 
momentum of the photon is subtracted from the momen¬ 
tum of the electron, i.e. the electron is recoiled, imposing 
the conservation law Xe ^ Xe ~ X 7 EO]- The simula¬ 
tion then proceeds by propagating the electron (via the 
Lorentz equation) and the photon (on a linear trajectory) 
to the next time step. In this way, multiple emissions 
are described as sequential single photon emissions, as in 
([7]), occurring at discrete time intervals (for further in¬ 
formation on the multi-photon calculation see Ref. [3J| 
and for discussions of double emissions Refs. HH E]). 
This method has recently been tested against cases where 
we can calculate the Compton spectra analytically and 
found to perform extremely well [42] , 


IV. RESULTS 
A. Plane wave model 

We begin by considering an electron with an initial 
70 = 800 in a head on collision with a plane wave laser 
of peak intensity ao = 200, A = 0 . 8 /im and duration 30fs 
FWHM. The emission spectra, calculated both classically 
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FIG. 4. Density plot showing how the electron 7 -factor 
changes with time (statistical distribution generated by 
recording the paths of 500 QED electrons all with the same 
initial condition 70 = 800). Parameters are ao = 200, 
A = 0 .8/im and duration 30fs. The white line shows the 7 - 
factor for a classical electron. 

and using the statistical QED routines, are plotted in 
Fig-1 It can be seen that, although for these parameters 
there is little difference between the frequency spectra, 
the classical theory predicts a much stronger signal at 
small angles than the QED theory. 

To understand the reason for this discrepancy we must 
consider how the electron energy changes with time. In 
Fig. [4] we plot the 7-factor of the classical electron (white 
line). It can be seen that the particle rapidly loses en¬ 
ergy as it approaches the pulse focus. This means that 
the angular direction of the emitted radiation will be con¬ 
tinuously changing as the ratio of the electron 7-factor to 
laser ao changes with time (see Fig. [TJ. This broadening 
of the angular range was first proposed in Ref. [43] as 
a signature of classical RR effects. (We also note that 
the rapid energy loss due to RR results in a natural re¬ 
sistance to the high-energy high-intensity regime, and is 
what prevents the use of lasers for distinguishing between 
RR models, see e.g. [44].) Superimposed on the same 
plot we also show how the electron 7-factor changes in 
the QED case. To generate this statistical density we 
ran the code 500 times, using the same initial conditions 
for each run. It can be seen clearly that the classical 
expressions overestimate the energy loss of the particle. 
The reason for this can be understood from the fact that 
the QED electron only emits at discrete times, allowing 
it to penetrate deeper into the laser pulse before it loses 
energy from an emission (an effect known as “straggling” 
[45]). whereas the classical electron is radiating and losing 
energy continuously [28, 39]. This means that, as the 


2 We note also that in Refs. |2llli6l it has been shown formally 


FIG. 5. Details of the setup for a more realistic example. 
Top left panel: We consider a bunch of 1000 electrons ran¬ 
domly distributed in transverse space according to a Gaussian 
distribution of 5jum FWHM. Top right panel: initial energy 
distribution of the electron bunch. Bottom panel: plot show¬ 
ing the laser intensity (colorscale shows ao) together with the 
electron bunch in the z — x plane. The solid black lines show 
how the laser waist varies as the pulse propagates (note that 
in the simulations the collision occurs at the origin). 

classical particle propagates through the pulse, its energy 
is going to be correspondingly lower than its QED coun¬ 
terpart. From Fig. [l] we can see that this means that the 
classical emissions will be at a correspondingly smaller 
angle than in the QED case. This is why the QED emis¬ 
sion spectrum dies off at a bigger angle than the classical 
spectrum. 

B. Realistic Example 

To show that these results still hold in an actual ex¬ 
periment we consider a more realistic setup. Instead of a 
plane wave field we now model our laser pulse as a 30fs 
duration paraxial beam, of wavelength A = 0.8/im, peak 
ao = 250 (equivalent to 2.64 x 10 23 W/cm 2 ), focussed to 
a waist radius of 2.5/im. These fields provide an accurate 
description of a focussed laser pulse, satisfying Maxwell’s 
equations to the order of the paraxial expansion param¬ 
eter #0 = ^/ttwo ~ 0.102, for more details see, e.g., 
Refs. [24j|47] (and for a discussion of Compton scattering 
in ultra-short focussed pulses see Ref. [48]). The parame¬ 
ters we have chosen are typical of what will be achievable 
at the ELI facility [5j. We model our electron source as 


that, for any given frequency, the classical radiation spectrum 
forms an upper bound to the (single photon) QED Compton 
spectrum. 
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FIG. 6. Emission spectrum for the realistic case (ao = 250) 
described in Fig. [5] The centre panel shows the radiation 
intensity as a function of frequency and angle. This panel is 
split into two, the top half showing the emissions for the QED 
simulation and the bottom half the classical. The right hand 
panel shows the total angular rate summed over all frequencies 
(both classical and QED for all angles), and the bottom panel 
the total frequency rate summed over all angles. Red lines: 
QED. Blue lines: classical. 


large ao in order that RR effects are significant and, ad¬ 
ditionally, one needs for the electrons to become reflected 
roughly in the middle of the pulse. If one has ao 7o 
then the electrons will be reflected very early and thus 
spend a significant amount of time co-propagating with 
the pulse - a regime in which the dynamics will be wholly 
classical. Also, if one has 70 ao the electrons will not 
be reflected and the angular broadening due to RR will 
be small. Nevertheless, there is a broad range of param¬ 
eters where ao and 7 are of same order of magnitude and 
the effect will be significant. (This parameter regime is 
found to be optimal for various studies of intense laser- 
particle interactions, including nonlinear Compton scat¬ 
tering ED, classical radiation reaction [43] , and for the 
generation of attosecond gamma-ray pulses [49].) 

The angular narrowing we have predicted is impor¬ 
tant for two reasons. Firstly, it provides a clear signal of 
strong field QED effects, distinct from classical RR, at 
parameters that will soon be obtainable. Secondly, it is 
an effect which will have to be taken into consideration 
when planning applications of Compton scattering with 
the new generation of ultra-intense lasers. 
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a beam of particles initially following a Gaussian distri¬ 
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The resulting emission spectra of the bunch of 1000 
electrons are shown in Fig. [6} Even though we are now us¬ 
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clearly evident. The classical theory predicts a strong 
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these signatures are sufficiently robust to survive realistic 
beam focussing effects. 


V. DISCUSSION 


This research was supported by the Swedish Research 
Council, grants # 2012-5644 and 2013-4248 and by the 
Wallenberg Foundation grant “ Plasma based compact ion 
accelerators". A.G. also acknowledges the Russian Foun¬ 
dation for Basic Research, project No. 15-37-21015. 


Appendix: Deriving the classical rate expression 
from that used in the QED case 

Here we demonstrate explicitly that the expression for 
photon emission in the QED form J7I coincides exactly 
with the one in the classical form (EF in the limit of low 
intensity and/or energy. Following the derivation of the 
classical synchrotron spectrum m we represent Eq. Q 
via the first and the second synchrotron function [50[ . 
Firstly, we convert the rate of emission to the spectral 
intensity of emission, multiplying it by photon energy a/ 
(as previously, according to the choice of units h = c— 1) 


Our analysis shows a significant difference in the angu¬ 
lar emission spectra predicted by the classical and QED 
theories. The intensities we have chosen are only slightly 
higher than the current state-of-the-art, and well within 
the parameters expected at the new generation of facil¬ 
ities, such as ELL For the regimes we have considered 
we find % e ^ 0.2 — 0.25, indicating that this is an ef¬ 
fect which becomes important before the onset of more 
explicit QED processes, such as pair production and run¬ 
away cascading. The difference in angular spectra won’t 
be visible for all sets of parameter values. One needs a 


dl q dT J 


duj' dX'yduj/dX'y 


A/37T7 2 


Y^) K VS (x)-/ dyK 1/3 (v) 

(A.l) 


Next, we apply a recursive relation of the modified Bessel 
function 


2 (t) = -K v -1 (t) - K v+ 1 0 1) (A.2) 
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for v = 2/3 to transform the integral in the right side 
(note, Kp(t) = K-n(t)) 


Using the ratio of the photon energy to the initial electron 
energy d = a// (my) = x 7 /y e as a parameter we obtain 


dl q uo'e 


/ „2 


<9o/ v^Try 2 


l+x 


^2/3 (A) 


+ J dy (^K 5 /3 (y) + 2—K 2 /3 (yfj ■ 


(A.3) 


Taking into account that lim^oo K^/^fy) = 0 we obtain 


dl q _ VZuj’e 2 Xe (Xe - X 7 ) 


dcu 1 


2irx 2 


X-y 


l + x 


F2 (x) + -Fi (x) 


(A.4) 

where F 2 (£) = £^ 2 / 3 (0 is the second synchrotron func¬ 
tion. In terms of the effective magnetic field H e ff we can 
express 


Xe = 7 


#eff 

E r , ’ 


x 7 = 


m E cr 


(A.5) 


<9o/ 


\/3 e s H e ff 
27r m 


(1 - d) x 


JL F 

l- d 2 V3xe 1 - dy 


+ Fi 


2 d 
3Xe 1 - d 


(A.6) 


This form of the expression makes it possible to see that 
in the case where the photon has much lower energy than 
the electron (d <C 1) the expression tends exactly to the 
classical one 
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